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We interpret the observed radial- velocity curve of the optical star in the low-mass X-ray binary 
2S 0921-630 using a Roche model, taking into account the X-ray heating of the optical star and 
screening of X-rays coming from the relativistic object by the accretion disk. Consequences 
of possible anisotropy of the X-ray radiation are considered. We obtain relations between the 
Q\ | masses of the optical and compact (X-ray) components, m v and m x , for orbital inclinations 

i = 60°, 75° 90°. Including X-ray heating enabled us to reduce the compact object's 
^j. ■ mass by ~ 0.5 — 1M , compared to the case with no heating. Based on the K0III spectral 

type of the optical component (with a probable mass of m v ~ 2.9M ), we concluded that 
m x ~ 2.45 - 2.55M (for i = 75° - 90°). If the K0III star has lost a substantial part of its 
mass as a result of mass exchange, as in the V404 Cyg and GRS 1905+105 systems, and its 
mass is m v ~ 0.65 — O.75M , the compact object's mass is close to the standard mass of a 
neutron star, m x ~ 1.4M (for i = 75° — 90°). Thus, it is probable that the X-ray source in 
^ ■ the 2S 0921-630 binary is an accreting neutron star. 



INTRODUCTION 



Currently, mass estimates are available for ~ 30 neutron stars. The most accurate 
of these are the masses of radio pulsars in Hulse- Taylor binaries [1] and of the radio pulsar 
J0737-3039 [2]. According to these estimates, the masses of neutron stars are confined to 
a fairly narrow range, 1.25 — 1.44M . Nevertheless, theoretical arguments suggesting the 
possible existence of massive neutron stars with m^s — 2 — 3M Q have accumulated. 

First, theoretical computations indicate that the expected mass of a neutron star 
formed during the core collapse of a massive star can lie in the range 1.0 — 1.8M [3]- [5]. 

Second, there exist quite a few hard equations of state for neutron matter for which 
the Oppenheimer-Volkoff mass exceeds 1.8M Q [6]. We especially note in this context 
recently published papers on so-called "Skyrmion stars". In 1999, Ouyed and Butler [7j 
considered an equation of state based on the model of Skyrme [8]. A characteristic feature 
of neutron-star models based on Skyrme's equation is their high upper mass limit: 2.95M 
for nonrotating and 3.45M for rotating stars 0, [TO] . 

Third, evolutionary computations of neutron-star masses in binary systems using the 
"Scenario Machine" [UJ have demonstrated that there could exist binary-system evolu- 
tionary channels, for which a neutron star is able to increase its mass via accretion by 
more than ~ 1M [T2l H5]. 

Because of these arguments, special interest has been attracted by the compact object 
in the low-mass X-ray binary 2S 0921-630 This binary X-ray system consists of a 
compact object and a low-mass optical star (V395 Car) [15] of spectral type KOIII [16] . 
The binary's orbital period is P or h = 9 d .006 ± 0.007 [17J. The nature of the compact 
object remains unknown: the binary does not exhibit phenomena characteristic of an 
X-ray pulsar or type 1 X-ray burster. Partial optical and X-ray eclipses are observed for 
the system, testifying to a high orbital inclination: i ~ 70° — 90° [16], [18] . 

According to Shahbaz et al.[16j, if the orbital inclination is i = 70° — 90°, the compact 
object's mass is 2.0 — 4.3M . The component-mass ratio, derived from the rotational line 
broadening, is q = m x /m v = 1.12 ± 0.18 [IB] . 

Jonker et al. [IT] present the results of fitting an accurate radial-velocity curve. As- 
suming i = 60° — 90°, the mass of the compact object was found to be 1.90 ± 0.25 < 
m x < 2.9 ± O.4M [17] . Given strong X-ray heating, these authors applied a so-called 
K-correction [T9] to the semiamplitude of the radial-velocity curve in order to estimate 
the component-mass ratio, q = m x /m v = 0.75 ± 0.37 [17] . 

The strength of the heating of the optical stars atmosphere by X-rays from the compact 
object k x = L x /L v ~ 10 complicates the fitting of the radial- velocity curve for the 2S 0921- 
630 system. The model of the binary system used to analyze the observed radial- velocity 
curve must take into account a variety of physical phenomena related to the compact 
objects high X-ray luminosity. 

Keeping in mind the importance of reliable mass estimates for the compact object 
in the 2S 0921-630 binary, we analyzed the high-precision data of [17] using the Roche 
model of [201 12] > taking into account the X-ray heating of the optical star, the screening 
of X-rays by the accretion disk, and the possible anisotropy of the X-ray radiation from 
the accretion disk associated with a black hole. 
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OBSERVATIONAL MATERIAL 



We used the spectroscopic data of jTTj for our analysis. These data were acquired 
between December 2003 and March 2004 with the VLT (Very Large Telescope). A total 
of 44 spectra with exposure times of 1300 s were obtained. The width of the spectrograph's 
slit was .4, making it possible to take spectra with a high resolution (0.75 A per pixel). 
The calibration uncertainty for the wavelength scale was 0.03 A. 

Spectroscopic data were obtained using two gratings, 1200R+93 and 1028z+29, which 
covered 5920-6520 A and 8360-8900 A respectively. The analysis of the spectrograms of 
the optical star in the 2S 0921-630 system performed in [17] confirmed its spectral type 
to be K0III. Radial velocities were determined by cross-correlating the spectra taken with 
the 1200R+93 diffraction grating with the spectra of a standard star. Stars with spectral 
types between G5 and K7, whose spectra were obtained with the Keck telescope with the 
same resolution, were used as standard stars. We adopted the 

We adopted the middle of the eclipse of the X-ray component by the optical star 
to be zero phase when plotting the radial-velocity curve. Jonker et al. [17] used their 
accurate radial-velocity curve to refine the binary's orbital period, which was found to 
be P orh = 9 d .006 ± 0.007 relative to the zero ephemeris JD Q = 2453000.49. The system's 
spatial velocity is 7 = 44.4 ±2.4 km/s. The semiamplitude of the observed radial- velocity 
curve, obtained using the K1V radial- velocity standard star HD 124106, is K v = 99.1 ±3.1 
km/s. The observed radial- velocity curve relative to the 7- velocity is shown in Fig{TJ 

ANALYSIS OF THE RADIAL- VELOCITY CURVES 

The optical star in the 2S 0921-630 close binary fills its Roche lobe. The star's 
surface is tidally and rotationally distorted, and the side facing the relativistic component 
is heated by X-rays from this component. The effects of the component's interaction must 
be taken into account when analyzing the optical component's radial-velocity curve. We 
fit the observed radial-velocity curve using the Roche model, taking into account the X- 
ray heating of the optical star. The algorithm we used is described in detail by Antokhina 
et al. [20] - [23]. We briefly summarize here the basic features of this method. 

The binary consists of an optical star treated using the Roche model and a point 
X-ray source. The optical star's surface is subdivided into several thousand (in our case, 
~ 2600) area elements, each of which we compute the emergent local flux assuming 
LTE, including the effect of the incident X-ray flux. Each area element is described by 
local temperature 7] oc , local gravity gi oc , and a parameter k l ° c , equal to the ratio of the 
incident X-ray flux to the emergent radiation flux, without taking into account external 
irradiation of the atmosphere. For these parameters at a given point on the surface, we 
solved the radiative-transport equations in a spectral line in the presence of the incident 
X-ray radiation, in order to compute a model atmosphere and obtain the intensity of the 
emergent radiation in the line and continuum. Our model synthesis of the radial-velocity 
curve was performed for the Cal 6439.075 A absorption line (hereafter, the Cal 6439 A 
line). For a given phase of the orbital period, we summed the contributions of the areas to 
the total flux, allowing for Doppler effects and the visibility of the areas for the observer. 
In this way, we computed the integrated continuum radiation flux from the star towards 
the observer, as well as the rotation-broadened profile of the spectral line, which was used 



2 



Table 1: Numerical parameters used to model the optical component's radial-velocity 
curve in the Roche model. 



P(day) 


9.006 


Orbital period 


m v (M Q ) 


var* 


Optical star's mass 


e 


0.0 


Eccentricity 


<) 


60, 75, 90 


Orbital inclination 


n 


1.0 


T-» 111 Oil* f* J f J 1 J * 1 J 

Roche lobe falling factor tor the optical component 


f 


1.0 


Rotational asynchronism factor for the optical component 


-Teff(K) 


4700 


Optical component s ettective temperature 


p 


0.08 


Gravity-darkening coefficient 


k x 


var* 


Ratio of the relativistic component's X-ray luminosity to the 
optical components bolometric luminosity, L x /L v 


A 


1.0 


Coefficient for reprocessing of the X-ray radiation 


u 


0.3 


Limb-darkening coefficient 


a p 


1.2 


Photon index of the X-ray spectrum 



* - These parameters of the X-ray binary were varied in our model fitting. 



to derive the star's radial velocity. The radial velocity of the star at a given orbital phase 
was determined from the mean wavelengths at residual intensity levels for the integrated 
absorption-line profile of one-third, one-half, and two-thirds. 

We carried out our analysis both including and excluding the effect of the instrumental 
profile on the model integrated profile of the Cal 6439 A line. The full width at half- 
maximum of the instrumental profile was taken to be FWHM=0.5A. The two sets of 
results were virtually identical, and we accordingly present here only those including the 
effect of the instrumental profile on the integrated profile of the Cal 6439 A line. 

We noted above that considerable X-ray heating of the optical star's atmosphere was 
observed in the binary system. Thus, the X-ray spectrum must be taken into account 
as accurately as possible, using the known X-ray spectrum [21]. Based on observations 
obtained with the XMM and Chandra observatories [24] . we adopted a photon index for 
the power-law spectrum a p = 1.2 in the range 0.1-12 keV. 

The compact object's intrinsic X-ray luminosity indicated by the XMM and Chandra 
data is L x ~ 10 36 erg/ s [24J . Assuming that the mean effective temperature of the optical 
star filling its Roche lobe is T e ff = 4700 K, its luminosity is L v ~ 2 ■ 10 35 erg/s. Given 
the uncertainties in the estimates of the distance to the system and of L x , we can assume 
k x = L x /L v ~ 10. We fit the observed radial- velocity curve using the atmospheric heating 
coefficients k x = and k x = 10. 

We fit the observed radial-velocity curve in three different models: with screening of 
the X-rays by the accretion disk both taken and not taken into account, and by allowing 
for the anisotropy of the X-rays emitted by the accretion disk of a rotating black hole. 
The numerical parameters of the Roche model for the 2S 0921-630 binary are presented 
in Tabled] 
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Table 2: Dependence of mx on mv for the Roche model including isotropic X-ray heating 
of the optical component with k x = 10 and i = 60°, 75°, 90°. (without screening of the 
X-rays by the accretion disk). 



m v (M Q ) 


i = 60° 
m x (M Q ) 


i = 75° 
m x {M & ) 


i = 90° 
m x {M Q ) 


1.0 


2.00 


1.70 


1.55 


1.9 


2.55 


2.15 


2.00 


2.4 


2.80 


2.35 


2.25 


2.9 


3.00 


2.55 


2.45 



Model 1. Analysis of the Observed Radial-Velocity Curve without Screening of the X-ray 
Flux by the Accretion Disk 

Since the exact mass of the optical star is not known, we treated the masses of both 
binary components as parameters to be determined. Since the inverse problem is rather 
cumbersome in this formulation, we carried out an exhaustive search for the parameters, 
thereby repeatedly solving the direct problem. For each mass of the optical component m v 
from the discrete set of values 1.0, 1.9, 2.4, 2.9M Q , fixing i, we carried out an exhaustive 
search for the compact object's mass m x . We tested the model's fit to the observations 
using the usual \ 2 statistical test, working with the significance levels a = 1% and 5% 
(cf. [25J for details). 

To quantitatively estimate the effect of X-ray heating on our estimates of m x , we fit the 
observed radial- velocity curve with X-ray heating both taken and not taken into account, 
and also allow for possible anisotropy of X-ray radiation from inner parts of the accretion 
disk around a black hole. When X-ray heating was taken into account, we set k x = 10. 
We obtained fits for inclinations i = 60°, 75° 90°. This analysis yields relations between 
the masses of the optical and compact objects, m v and m x . The numerical results are 
collected in Tables [2] and [3] and are plotted in Fig. [2j 

Tables [2] and [3] present the values of m x corresponding to the minimum x 2 residuals for 
the difference between the theoretical and observed radial-velocity curves. The values of 
m x in Tables [2] and 33] are not accompanied by uncertainties, because both of the binary 
models were rejected at the a = 1% and 5% significance levels. The quantiles for these 
significance levels are A x % = 38.93 and A 5 % = 32.67. When X-ray heating of the optical 
star is included (k x = 10), the minimum residual is x m in — 120 — 150. When it is not 
included (k x = 0), the minimum residual increases to Xmin - 350 - 380 - The statistical 
inadequacy of the models in describing the observational data is due to the considerable 
scatter of the data points, with uncertainties in the radial velocities for each data point 
being comparatively small (see Fig. H] below). As was noted in [17], the considerable 
scatter of the data points in the radial-velocity curve is probably due to variable X-ray 
heating of the star. 

Let us compare the results of our fitting of the observed radial- velocity curve obtained 
using the Roche model both with (Table [2]) and without (Table [3]) X-ray heating. With 
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Table 3: Dependence of mx on mv for the Roche model without heating the optical 
component (k x = 0) for i = 60°, 75°, 90°. 



mJMn) 


i = 60° 
m x (M Q ) 


i = 75° 
m x (M Q ) 


i = 90° 
m x {M Q ) 


1.0 


2.55 


2.10 


1.95 


1.9 


3.30 


2.75 


2.60 


2.4 


3.65 


3.05 


2.85 


2.9 


4.00 


3.35 


3.15 



X-ray heating {k x = 10), with i = 60° — 90° and m v = 1.0 — 2.9M , the mass of the 
relativistic object is m x = 1.55 — 3.0M Q ; without X-ray heating (k x = 0), the compact 
object's mass is m x = 1.95 — 4.OM . Thus, the fit of the data without X-ray heating of 
the optical star's atmosphere leads to masses of the compact object that are too high by 
0.5 - l.OM 

Let us consider in more detail why m x is systematically too high in the model without 
X-ray heating. Figure [3] displays the model integrated profiles of the Cal 6439 A line at 
orbital phase <fi = 0.25 obtained for m v = 2AM & , m x = 2.35M , and i = 75°, with X-ray 
heating both taken (solid curve) and not taken (dashed curve) into account. 

When X-ray heating is taken into account, an emission component appears in the 
integrated profile of the absorption line (Fig. [3]). Because of this emission component, 
the "center of gravity" of the absorption line is displaced. Thus, when the X-ray heating 
is included, the line's center of gravity in Fig. [3] is at a wavelength of 6441.131 A, corre- 
sponding to a radial velocity of 95.44 km/s. If there is no X-ray heating, the line's center 
of gravity is at 6440.818 A, corresponding to a radial velocity of 81.57 km/s. Thus, due 
to the distortion of the absorption line profile by the emission component, the observed 
amplitude of the optical component's radial- velocity curve can be explained using a lower 
mass for the compact object. For this reason, carefully taking into account the X-ray heat- 
ing in the 2S 0921-63 binary reduces the compact object's mass from m x = 1.95 — 4.OM 
to m x = 1.55 - 3.OM . 

This is illustrated by model radial- velocity curves in|U For X-ray heating with k x = 10, 
m v = 2.4M , % = 75°, the minimum residual (Xmin = 149.8) is achieved for a mass of the 
compact object m x = 2.35M (solid curve in Figf4]). The semiamplitude of the optical 
componentfs radialvelocity curve is K v = 96.81 km/s. For the same binary parameters 
but without X-ray heating, the radial- velocity curve has a semiamplitude of K v = 81.27 
km/s (dash-dotted curve in Figf4]). For the analysis without X-ray heating and with 
m v = 2.4M and % = 75°, the minimum residual (Xmin = 356.7) is achieved for the higher 
mass of the compact object m v = 3.O5M (dashed curve in FigH]). 

Let us also consider the shape of the optical star's absorption line in the presence 
of strong X-ray heating. It is obvious that estimates of the mass of the compact object 
depend directly on the uncertainty in the center of gravity of the absorption line in the 
optical star's spectrum. When analyzing the shape of model integrated profiles for the 
Cal 6439 A line, we encountered the problem of determining the line's center of gravity. 
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Table 4: Dependence of m x on m v for the Roche model with i = 75°, heating of the 
optical component with k x — 10, and screening of the X-rays by the accretion disk. 



m v (M Q ) 


m x (M & ) 


1.0 


1.70 


1.9 


2.15 


2.4 


2.40 


2.9 


2.60 



For example, in a binary with m v = 2.4M , m v = 2.35M , and i = 75°, the center 
of gravity of the line at phases 0.47-0.50 cannot be determined unambiguously. Figj5] 
presents integrated profiles of the Cal 6439 A line for orbital phases 0.45-0.50. We can 
see that, beginning with orbital phase 0.47, the line profile has two minima, preventing 
unambiguous determination of the center of gravity of the line, and hence of the radial 
velocity. For this reason, a break in the radial-velocity curve is observed at phases 0.47 
and 0.53 (Fig. Sj). 

Because of the ambiguity in the absorption-line profiles in X-ray binary systems whose 
optical stars are heated, with k x > 10, it is difficult to derive a single satisfactory radial- 
velocity curve. Consequently, the masses of the relativistic objects must be estimated 
directly from the orbital variations of the observed line profiles, rather than from indirect 
data, such as radial- velocity curves. Modern 8- 10m telescopes can provide the needed 
high quality observational data for this purpose. 

Model 2. Analysis of the Observed Radial- Velocity Curve with Screening of the X-ray Flux 
by the Accretion Disk 

In this section, we allow for screening of the X-ray flux by the accretion disk, assuming 
the total screening angle of the accretion disk to be 5° and the accretion disk to be 
situated in the binary's orbital plane. Some area elements on the optical star happen to 
be in the strip that is shielded from the X-rays from the relativistic object. We computed 
absorption-line intensities for these areas proceeding as above, taking the local X-ray 
heating coefficient to be k l ° c = 0. We analyzed the observed radial- velocity curve for 
discretemasses of the optical star m v = 1.0, 1.9, 2.4, 2.9M , and for i = 75°. The results 
are presented in Table HI 

The results were quite close to (for masses m v = 1.0), and 1.9M , identical with those 
obtained in the Roche model without screening the X-ray flux by the accretion disk (Table 
[2]). For this reason, we did not continue our analysis for orbital inclinations i = 60° and 
90°. The shadow strip on the optical component turned out to be too narrow to have any 
considerable influence on the results. 

Similarly, we analyzed the observed radial- velocity curve allowing for screening of the 
X-ray flux by the accretion disk for the case k x = 30, with i = 75°. and m v = 1.0, 1.9, 
2.4, 2.9M , as in the previous case. The results are presented in Table [5j The derived 
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Table 5: Dependence of m x on m v for the Roche model with k x = 30 and i = 75°. 



m v (M Q ) 


with screening of the 
X-rays m x (M Q ) 


without screening of the 
X-rays m x (M & ) 


1.0 


1.65 


1.60 


1.9 


2.08 


2.05 


2.4 


2.30 


2.28 


2.9 


2.50 


2.45 



masses are tabulated without error intervals, because all the binary models were rejected 
at the 1% and 5% significance levels. The m x values given correspond to the minimum 
X 2 residuals. 

It follows from Table |5] that including screening of the X-ray flux L x = 30L V by the 
accretion disk has little influence on the results. The largest mass difference revealed for 
the two models is O.O5M (Table [5]). The effect of screening will be even smaller if the 
heating coefficient is lower. Thus, the X-ray pulsar masses we derived earlier in [27J that 
did not allow for screening of the X-ray flux by the accretion disks can be considered 
reliable. 

Analysis of the Observed Radial- Velocity Curve Allowing for Anisotropic X-ray Radiation 
from the Accretion Disk 

To complete our study, we considered the hypothesis that the 2S 0921-630 binary 
contains a low-mass black hole. In this case, the X-rays from the accretion disk should be 
anisotropic, reducing the X-ray flux in the orbital plane. Note that the presence of strong 
emission components for many absorption lines in the optical component's spectrum [T7] 
suggests that the X-ray radiation is isotropic. The emission lines in the spectrum of the 
optical star, which indicate considerable X-ray heating of the optical component, also 
provide indirect evidence that the system contains a neutron star. 

We assumed the accretion disk to be optically thin and located in the binary's orbital 
plane. The anisotropy of the X-rays from the accretion disk was described with Eqs. ([T]) 
and (ED from [28]: 

dL L x - F(9) 

dn 4vr ' 1 j 

F(9) = -cos 0(1 + 2cosfl), (2) 

where 9 is the angle between the normal to the disk plane and the direction of an element 
of solid angle dfl. The geometry of the binary system for the Roche model with the 
parameters in Table Q] suggests that the X-ray flux falls on the optical component at 
9 ~ 70° — 90°. We can see from fl2]) that the X-ray flux is considerably attenuated for 
such values of 9. Figure [6] presents the model integrated profiles of the Cal 6439 A line 
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Table 6: Dependence of m x on m v for the Roche model with k x = 10, anisotropic radiation 
from the X-ray source, and i = 60°, 75°, 90°. 





i = 60° 
m x (M Q ) 


i = 75° 
m x (M Q ) 


i = 90° 
m x (M & ) 


1.0 


2.20 


1.80 


1.70 


1.9 


2.75 


2.30 


2.15 


2.4 


3.00 


2.50 


2.35 


2.9 


3.25 


2.75 


2.60 



at orbital phase 0.25. Note that the emission component of the integrated line profile is 
weaker than in the case of an isotropic X-ray flux. 

The results of our fitting of the observed radialvelocity curve for the case of anisotropic 
X-rays from the accretion disk are presented in Table [6] and plotted in Fig. [71 We give 
the central values without uncertainties, since the binary models were rejected at the 
a = 1% and 5% significance levels. For m v = 1.0 — 2.9M Q , the mass of the black hole 
is 1.70 — 3.25M . Recall that the compact object's mass for the case of isotropic X-rays 
with k x = 10 was m x = 1.55 — 3.OM . Thus, allowing for an anisotropic X-ray flux from 
the accretion disk increases the compact object's mass by ~ O.2M compared to the case 
of an isotropic X-ray flux. 

DISCUSSION 

The results of our analysis of the high-precision radial- velocity curve p2] (Fig. [TJ) with 
orbital inclinations i = 60° — 90° and the optical component's mass m v = 1.0 — 2.9M Q 
are presented in Table [71 We consider the mass of the compact object obtained in Model 
1 to be our main result. 

As we noted above, the binary exhibits X-ray eclipses, testifying to high orbital in- 
clination: i ~ 70° — 90° [IE]. Since there are no X-ray dips, we can assume i > 80° 
|18j . For our mass estimates, we adopted values i = 75° — 90°. In this case, assuming 
m v = 1.0 — 2.9M , the mass of the compact object is m x ~ 1.55 — 2.55M (Fig. [2]). 
Taking into account the optical component's spectral type, K0III {m v ~ 2.9M [29J ) , we 
find the compact object's mass to be m x = 2.45 — 2.55M (Fig. [2]). 

The actual mass of the optical star in an interacting binary, m v , can differ from the 
estimate obtained based on its spectral type. For example, the mass estimates for the 
low-mass transient X-ray binaries V404 Cyg (spectral type of the optical star K0IV) and 
GRS 1915+105 (spectral type of the optical star Kill) derived from rotational broadening 
of absorption lines are O.7±O.1M |30j and 0.81 ±O.53M (21], respectively. At the same 
time, mass estimates based on the corresponding spectral types and luminosity classes 
are ~ 1.3M and ~ 2.3 — 2.9M [29], respectively. Thus, the estimates based on spectral 
type are too high by almost a factor of ~ 2 compared to the spectroscopic m v values. 
This is in agreement with current theoretical ideas about stellar evolution in low-mass 
close binaries. According to computations using the "Scenario Machine" [UJ, an optical 
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Table 7: Mass of the compact object in the 2S 0921-630 X-ray binary for various models, 
with k T = 10. 





m x (M & ) 


1 




1.55-3.0 


2 




1.55-3.0 


3 




1.70-3.25 


* - Model 1: isotropic X-ray 


source, 


without screening of 



the X-ray flux by the accretion disk; model 2: isotropic X- 
ray source, with screening of the X-ray flux by the accretion 
disk; model 3: anisotropic X-ray source, without screening 
of the X-ray flux by the accretion disk. 

component with initial mass ~ 3M loses more than ~ 1M Q during the mass-exchange 
stage [T2] . 

The optical stars in the X-ray binaries V404 Cyg, GRS 1915+105, and 2S 0921-630 
have similar spectral types and luminosity classes [T6"| l3~Tj . Assuming that the K0III star 
in the 2S 0921-630 system lost a substantial fraction of its mass via exchange and that 
its mass, like the optical star in GRS 1905+105, is m v ~ 1M , we find the mass of the 
compact object in the 2S 0921-630 system to be m x ~ 1.6 - 1.7M (for i = 75° - 90°). 
This is close to the mean mass of a neutron star, making it likely that we are dealing with 
an accreting neutron star in the 2S 0921-630 binary. 

To extend our analysis for the 2S 0921-630 binary, we considered the possibility that 
the compact object's mass is close to the standard mass of a neutron star, m x = 1AM Q . 
In this case, the mass of the optical component is m v = 0.65 — 0.75M Q , where we have 
assumed i = 75° — 90° since partial X-ray eclipses are observed for the system [TO] . 
Formally, from the stellar-evolution point of view, single stars with such masses cannot 
become giants over the Hubble time. In this case, the presence of a relativistic companion 
enables a star with an initial mass > O.8M to lose some its mass during the semidetached 
phase [32] . 

The result of the evolution of a giant with mass ~ 0.6 — O.7M is a helium white 
dwarf with a mass ~ 0.3 — 0.4M Q [32J. Thus, it is quite plausible that further evolution 
of the 2S 0921-630 X-ray binary will result in a millisecond radio pulsar in a binary with 
a helium white dwarf. This hypothesis is in good quantitative agreement with the model 
computations of |32j and the empirical relation between the orbital periods and secondary 
(white dwarf) masses for binaries with millisecond radio pulsars (cf. [[32], Fig. 4] and 
[[M], Fig. 9]). 

We note again that we rejected all models for the X-ray binary system at the a = 
1% and 5% significance levels. Tables [2] - [7] present the m x values corresponding to 
the minimum x 2 residuals. We must bear in mind that, in all cases, the models do 
not completely adequately represent the observational data, so that our estimates of m x 
cannot be considered final. 
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Note the importance of the methodological results of our study. Our computations for 
the Roche model demonstrated that the profile of the absorption line is rather complex 
in the case of strong X-ray heating, and suffers considerable variations in the course of 
orbital motion (Figs. [3J [5]). Thus, estimates of q based on rotational broadening of 
absorption-line profiles in low-mass X-ray binaries with strong X-ray heating cannot be 
considered trustworthy. Our future plans include a study of this problem. 

CONCLUSIONS 

Taking into account large orbital inclination, i = 75° — 90°, and our relations between 
the masses of the components in the 2S 0921-630 system (Fig. [2]), we conclude that 
the mass of the compact object is m x ~ 1.55 — 2.55M Q if the optical star's mass is 
m v = 1.0 — 2.9M . Our study demonstrates that, if we allow for the possible mass loss by 
the optical component that reduces its mass to m v = 1M , the presence of an accreting 
neutron star with m x ~ 1.6 — 1.7M is most likely for this binary. The observed X-ray 
spectrum [21] of the compact object in 2S 0921-630 could correspond to a black hole, 
as well as to a neutron star with a weak magnetic field. It remains difficult to make an 
unambiguous selection between these possibilities for the nature of the compact object. 

Our study demonstrates that high-resolution spectrograms (A/ A A = 50000 — 100000) 
are needed to estimate the compact object's mass. The most correct approach to estimat- 
ing this mass, and hence revealing the nature of the compact object, is to fit the orbital 
variations of absorption-line profiles in the spectrum of the optical component. 

We thank A.V. Tutukov for helpful comments. 
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Figure 1: Observed radial- velocity curve of the optical component of the X-ray binary 2S 
0921-630 from [T7j and the theoretical radial-velocity curve in the Roche model, for the 
compact-object mass m x = 2.35M , optical star mass m v = 2.4M , X-ray heating coef- 
ficient for the optical star k x = 10, and orbital inclination % = 75° (the remaining binary 
parameters are collected in Tabled])). The theoretical radial- velocity curve corresponds 
to the minimum-residual fit (Xmin = 149.8). 
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Figure 2: Relations between the masses of the optical star and the relativistic component in the 
2S 0921-63 binary from our fitting of the observed radial-velocity curve using the Roche model, 
assuming k x = 10 (solid broken lines) and k x = (dashed broken lines) for X-ray heating. The 
X-ray heating is taken to be isotropic. 
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Figure 3: Model integrated profiles of the Cal 6439 A line of the optical component in the 
2S 0921-630 X-ray binary at orbital phase 0.25. The profiles were obtained for the Roche 
model with m x = 2.35M , m v = 2.4M , and % = 75°. The dashed curve is the line profile 
derived in the Roche model without including the reflection effect. The solid curve is the 
line profile obtained with X-ray heating of the optical star corresponding to k x = 10. Both 
model integrated profiles were convolved with the spectrograph's instrumental profile, 
assuming its FWHM=0.5A. 
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Figure 4: Observed radial-velocity curve from [T7j and the theoretical radial-velocity 
curves obtained for the Roche model with i = 75. Shown are the best-fit radial- velocity 
curves for X-ray heating with k x = 10, with m v = 2AM & and m x = 2.35M Q (solid); 
without X-ray heating {k x = 0), with m v = 2.4M and m x = 2.35M (dash-dotted); and 
without X-ray heating {k x = 0), with m v = 2.4M and m x = 3.05M Q (dashed). The 
other parameters of the binary are given in Table [TJ 
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Figure 5: Model integrated profiles of the Cal 6439 A line obtained for the Roche model 
with X-ray heating of the optical component (k x = 10) and for m v = 2.4M , m x = 
2.35M , % = 75°. The other binary parameters are given in Table 1. Our notation for the 
orbital phases is shown near the curves. 
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Figure 6: Model integrated profiles of the Cal 6439A line of the optical component in 
the 2S 0921-630 X-ray binary system at orbital phase 0.25. The profiles were obtained 
for the Roche model with m x = 2.35M , m v = 2.4M , % = 90°, and k x = 10. Shown 
are the absorption-line profile assuming isotropic (solid) and anisotropic (dashed) X-ray 
radiation from the compact object. Both model integrated profiles were convolved with 
the spectrograph's instrumental profile, assuming its FWHM=0.5A. 
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Figure 7: Relation between m v and m x in the 2S 0921-630 binary obtained by fitting the 
observed radial-velocity curve for the Roche model with k x = 10 and anisotropic radiation from 
the accretion disk. See text for details. 



19 



